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The forced intrusion of water in hydrophobic nanoporous pulverulent material is of interest for 
quick storage of energy. With nanometric pores the energy storage capacity is controlled by inter¬ 
facial phenomena. With subnanometric pores, we demonstrate that a breakdown occurs with the 
emergence of molecular exclusion as a leading contribution. This bulk exclusion effect leads to an 
osmotic contribution to the pressure that can reach levels never previously sustained. We illustrate 
on various electrolytes and different microporous materials, that a simple osmotic pressure law ac¬ 
counts quantitatively for the enhancement of the intrusion and extrusion pressures governing the 
forced wetting and spontaneous drying of the nanopores. Using electrolyte solutions, energy storage 
and power capacities can be widely enhanced. 


Water in hydrophobic confinement is a topic of major 
importance for industrial applications, such as boiling, 
heat or mass transfers at interfaces nHa, but also for the 
general understanding of hydrophobic interactions medi¬ 
ated by water in biological matter ISHZI and biomolecular 
responses under osmotic stress [8]. Ordered hydropho¬ 
bic nanoporous materials have been used to study water 
confined between hydrophobic surfaces in well defined 
geometries. But these material are also considered for 
novel energy applications. Due to the large specific area 
of nanoporous materials, of the order of 1000 m^/cm“^, 
and to the high tension of liquid/hydrophobic solid in¬ 
terfaces, forcing the intrusion of water in these materials 
provides a way to store energy in the form of interfacial 
free energy [3] . Then, the spontaneous extrusion of water 
out of the nanopores, related to their drying transition, 
allows the partial or full recovery of the stored energy. 

In nanoporous materials of pore diameter larger than 
typically ten water molecules, it was shown that macro¬ 
scopic concepts describe quantitatively the pressure Pint 
at which the forced wetting, or intrusion, occurs. For 
instance in the cylindrical pores of hydrophobic meso- 
porous silicas, the intrusion pressure obeys the Laplace 
law of classical capillarity and scales as the inverse of 
the pore radius IjRp [9]. The drying transition on the 
other hand was shown to be triggered by the nucleation 
of nanobubbles, a mechanism also governed by interfa- 
cial phenomena pUUT3] . This nucleation process, in con¬ 
finement, has also proven to be of interest to study the 
growth of nanobubbles, of precisely controlled shape, in 
relation with line tension effects. 

With subnanometric pores, such as the pores of ze- 
olithes or Metal Organic Frameworks, a breakdown oc¬ 
curs m- In particular, electrolyte solutions exhibit huge 
wetting and drying pressures as compared to the ones of 
pure water. In this letter, we demonstrate for the first 
time that exclusion effects emerge as a leading contri¬ 
bution responsible for this giant pressure increase. We 


find that during the intrusion/extrusion processes, the 
pressure in the bulk electrolyte corresponds to the sim¬ 
ple addition of the osmotic pressure 11 to the value of 
the intrusion/extrusion pressures obtained for pure wa¬ 
ter. This simple osmotic effect can increase very signifi¬ 
cantly the density of the stored energy. 

Experiments have been carried out on a Metal Organic 
Framework (MOF) called Zeolitic Imidazolate Frame¬ 
work 8 (ZIF-8)[T5], combined with seven different elec¬ 
trolyte solutions in various concentrations and temper¬ 
ature. ZIF-8 is commonly studied in the field of gas 
separation because of its high selectivity [16]. An im¬ 
portant feature of ZIF-8 is its natural hydrophobicity. 
Forced intrusion/extrusion cycles were achieved using a 
piston-cylinder device as described by Guillemot et al 
[T7] (see figure [^, mounted on a traction/compression 
machine. Pressure-Volume (P-V) cycles were obtained 
on previously evacuated samples at constant tempera¬ 
ture and under quasi-static conditions. 





FIG. 1. Principle of forced intrusion/extrusion cycles. Pres¬ 
sure measurements are performed in an instrumented water¬ 
proof chamber, containing ZIF-8 (red) [15] and electrolytes so¬ 
lution (blue), mounted on a traction machine. 


The pressure of the liquid versus the intruded volume 
for a ZIF-8/water system is presented in the inset of fig¬ 
ure Intrusion and extrusion processes correspond to 
plateaus of the cycle. The value of the intrusion (resp. 
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extrusion) pressure is defined as the average value on the 
intrusion (resp. extrusion) plateau. The ZIF-8/pure wa¬ 
ter system shows a small hysteresis with flat plateaus 
for both intrusion and extrusion. The intruded volume 
is about 500 mm^.g“^. Values of the intrusion and ex- 
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FIG. 2. Values of the intrusion and extrusion pressures in 
ZIF-8 as a function of the concentration c of NaCl solutions 
at 323 K. The inset shows the difference between the pressure 
and the van ’t Hoff osmotic pressure versus intruded volume 
for ZIF-8/water and two different ZIF-8/NaCl solutions at 
T = 343 K. 


trusion pressures are respectively 25.2 ± 0.2 MPa and 
20.9 ±0.2 MPa for this pure water system. The cycle ob¬ 
tained for ZIF8/NaCl solutions reproduce perfectly the 
pure water cycle with a simple shift in pressure (see in¬ 
set of figure . As shown in figure the values of the 
intrusion pressure Pint and the extrusion pressure Pext 
increase linearly with the NaCl concentration c, with 
slopes of respectively 5.5 ±0.2 MPa.M“^ for intrusion and 
5.2 ±0.2 MPa.M“^ for extrusion. This linear dependency 
corresponds very closely to the van’t Hoff law of osmotic 
pressure H: 


H = icRT (1) 

where R is the perfect gas constant, T the temperature 
and i the number of ions per electrolyte molecule. The 
corresponding slope is iRT = 5.37 MPa.M“^ for i = 2 and 
T = 323K. This theoretical value is very close to the 
slopes of the experimental data. 

This behaviour unravels the mechanism at work in the 
intrusion-extrusion process of the NaCl solutions. The 
ZIF-8 behaves as an almost ideal sieve for Na”^ and CF 
ions, as predicted by molecular simulations [18]. Because 
of this selectivity, only water penetrates into the material, 
as sketched in figure The pressure difference between 


the pure water inside the material and the surrounding 
solution is the osmotic pressure. This pressure difference 
is constant through the whole process, therefore the elec¬ 
trolyte cycles are shifted replicas of the pure water cycles, 
as only pure water is present into the material. In other 
words, the pressure shift observed with the electrolyte so¬ 
lution is the additional pressure required to extract pure 
water from the surrounding electrolyte. 
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FIG. 3. Schematic view of the ion exclusion from the ZIF- 
8 microporous material leading to an osmotic contribution. 
The jump between the pressure pw inside the pores of a mi¬ 
croporous particle and the pressure ± H outside the pores 
is the osmotic pressure H. 


We expect the same scenario to hold for all cases in 
which ions are perfectly excluded from the microporous 
material. We have studied seven electrolytes: lithium 
chloride LiCl, lithium iodide Lil, sodium chloride NaCl, 
sodium bromide NaBr, sodium iodide Nal, calcium chlo¬ 
ride CaCl 2 and cesium chloride CsCl. For each solu¬ 
tion, the excess pressure with respect to pure water, 
AP = psoiution _ pwater determined both for in¬ 

trusion and for extrusion at different temperatures m- 
Figure summarizes the results with the experimental 
AP’s plotted versus the van’t Hoff pressure icRT. The 
figure also includes data obtained in Silicalite-l zeolite 
with NaCl and LiCl mi and data obtained in ZIF-8 with 
NaCl, LiCl and KCl [20]. The case of Nal and Lil is dis¬ 
cussed later and reported in figure As predicted by 
van’t Hoff law, data points obtained for low concentra¬ 
tions of LiCl, NaCl, NaBr and CsCl, align along a straight 
line whose slope is 1. 

Measurable deviations from the van ’t Hoff law are 
observed with saturating solutions of NaCl and LiCl in 
Silicalite-l zeolites (Tzanis data) as well as with CaCl 2 - 
0.5M in ZIF-8 (our data). The deviating pressures are 
higher than predicted by the van ’t Hoff law. We at¬ 
tribute these deviations to non-ideal effects in the elec¬ 
trolyte solutions. Non-ideal effects are expected when the 
Coulomb interaction energy between ions becomes com¬ 
parable to the thermal energy. This occurs for a typical 
distance between ions scaling as z^z~Xb with and 2 :“ 
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FIG. 4. Excess intrusion and extrusion pressures with re¬ 
spect to pure water measured for the electrolytes mentioned 
in the legend. The x-scale is the van ’t Hoff osmotic pres¬ 
sure n = icRT with i the number of ions per salt molecule 
and c the salt concentration. The values of our own data (■) 
are listed in supplementary Tables 1 to 5 m- {O) are data 
from Tzanis et al IH, (O) are data from Ortiz et al [20] . 
The straight line (•••) corresponds to the van ’t Hoff law of 
osmotic pressure. The continuous line (— ) is a simulation of 
the osmotic pressure of NaCl solutions from Luo et al [21] . 
The inset is a magnification of our data with (O) intrusion 
pressure and (A) extrusion pressure. 


the ions valence and Xb the Bjerrum length: 

~ ineolrksT 

with eo the vacuum permitivity, the water dielectric 
constant and ks the Boltzmann constant. The distance 
z'^z~Xb is related to a limit salt concentration Cz+z-'- 

_ 1 1 _ 1 /47reoerA:BT\^ 

and a limit pressure liz+z- = iCz+z-RT above which cor¬ 
relation effects between ions can not be neglected. For 
monovalent ions at the lowest considered temperature 
T = 303 K, the limit concentration is around cn = 1 M, 
corresponding to a limit pressure Hu of about 5.5 MPa. 
At concentrations much higher than cn the osmotic pres¬ 
sure is expected to increase more rapidly than the linear 
van’t Hoff prediction. This is the range of concentration 
above which deviations are experimentally observed for 
NaCl and LiCl in figure The molecular dynamic sim¬ 
ulation of the osmotic pressure of NaCl carried out by 
Luo et al m indeed meets the data point measured by 


Tzanis et al and Ortiz et al for saturating NaCl solutions 
in zeolithes [Til 120]. Taking into account the trend of the 
simulation, it is safe to propose that osmotic pressure is 
also the good candidate to explain the much higher pres¬ 
sure measured by the same authors for concentrated LiCl 
solutions (180 MPa of excess pressure at saturation). 

The valence of the ions has a major impact on liz+z - 
For a salt with one divalent ion and one monovalent ion 
such as CaCl 2 , the limit concentration is C 21 = 0.15 M, 
and the corresponding pressure is n 2 i = 1.2 MPa. Ac¬ 
cordingly, the experimental values of AP measured for 
CaCl 2 solutions are larger than the van’t Hoff theoreti¬ 
cal curve in figure and the data points go away from 
the van’t Hoff law as the concentration increases, in a 
similar way to what is observed with monovalent ions at 
higher concentration. 

Remains the particular case of Nal and Lil solutions 
depicted on figure The behaviour of these salts is very 
different from the previous ones. The cycle of the Nal - 
0.94 M solution is compared to one of pure water in the 
inset of figure The electrolyte cycle is not anymore a 
simple translation of the pure water cycle. The excess of 
the intrusion pressure is much larger than the excess of 
the extrusion pressure, and the cycle hysteresis is larger 
for the electrolyte solution. A similar result is obtained 
with all concentrations of Nal and Lil solutions as shown 
in figure Furthermore the excess pressures measured 
in the electrolyte cycle with respect to pure water are 
below the van’t Hoff osmotic pressure. 



FIG. 5. Excess intrusion pressure (O) and excess extrusion 
pressure (A) of Nal and Lil solutions in ZIE-8 with respect 
to pure water, plotted as a function of the van ’t Hoff osmotic 
pressure icRT. Inset : difference between the pressure and 
the van ’t Hoff osmotic pressure of ZIE-8/water (blue) and 
ZIE-8/NaI - 0.94M (red) at 343 K. 
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We attribute this behavior to the partial penetration 
of ions inside the micropores. As a result the interfacial 
energies are expected to be different from the one of pure 
water, and consequently the shape of the cycle is also dif¬ 
ferent. Furthermore the pressure difference between the 
surrounding electrolyte and the intruding liquid is less 
than the van ’t Hoff osmotic pressure as the intruding 
liquid is not pure water. The penetration of these ions 
into the pores is far from intuitive as the radius of the 
hydrated I“, vh = 3.3 A[22], is twice larger than the ra¬ 
dius of the smallest ZIF-8 pores. The partial penetration 
of 1“ ions in ZIF-8 is certainly due to their large polar- 
isability, responsible for their deformability and affinity 
for hydrophobic surfaces [23l [24] . 

In any case, the transition between ion intrusion and 
exclusion is not related to electrostatic screening. Iodide 
ions enter the pores for small concentrations correspond¬ 
ing to Debye length larger than the pore size while other 
ions are still excluded from the pore at high concentra¬ 
tions corresponding to Debye length smaller than 0.1 nm. 

In summary, our results show that electrolyte solutions 
can exhibit giant intrusion/extrusion pressures in micro¬ 
pores due to osmotic effects. As ions are repelled from 
the pores rather than confined in them, the intrusion pro¬ 
cess is a three-dimensional extraction of pure water from 
the electrolyte solution. This osmotic mechanism allows 
to increase the pressure up to one order of magnitude 
compared to the intrusion/extrusion pressure obtained 
with pure water. The intrusion process is associated to a 
storage of mechanical energy, equal to the product of the 
macroscopic pressure times the intruded volume [25]. As 
a result the energy stored in the pores per unit volume 
can be higher than with ultracapacitors. Morevover, it 
was demonstrated that the intrusion/extrusion pressures 
depend weakly on the rate of the process due the very 
short time scale of wetting and drying transitions [12]. 
Therefore our results pave the way towards the design of 
lyophobic energy storage systems with targeted power- 
densities as large as 50 kW per solution liter that is one 
order of magnitude larger than with ultracapacitors. 

From a more fundamental point of view, our results 
demonstrate experimentally that the concept of osmotic 
pressure can be extended up to giant values of hundreds 
of MPa that have never been reached with conventional 
separation membranes. The selective pulverulent nano¬ 
material act as a volumic membrane able to sustain huge 
osmotic pressures because each individual material par¬ 
ticle is submitted to a uniform osmotic pressure. Finally, 
the forced intrusion of aqueous solution in hydropho¬ 
bic nanoporous material can be put forward as a way 
to probe mechanically molecular interactions. This me¬ 
chanical characterization offers a new way to probe the 
stability of the hydration layer of ions and the affinity 
of ions for hydrophobic surfaces as exemplified with the 
difference of behaviors between iodide ions and other an¬ 
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